Dirac fermions in a power-law correlated random vector potential 



o 
o 

O 



CD 



o 

> 

in 
o 



in 
o 

o 



X 



D. V. Khveshchenko 

Department of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599 

We study localization properties of two-dimensional Dirac fermions subject to a power-law cor- 
related random vector potential describing, e.g., the effect of "ripples" in graphene. By using a 
variety of techniques (low-order perturbation theory, self-consistent Born approximation, replicas, 
and supersymmetry) we make a case for a possible complete localization of all the electronic states 
and compute the density of states. 



Anderson localization in two spatial dimensions (2D) 
has traditionally been studied in the context of doped 
semiconductors which are characterized by the presence 
of fermionic excitations with an extended Fermi sur- 
face. After having been maturing steadily for almost 
three decades, the theory of 2D localization was rekin- 
dled by the advent of high- Tc superconductors and, more 
recently, graphene, in which systems the Fermi surface 
of quasiparticle excitations consists of isolated (nodal) 
Fermi points. 

Unlike in the conventional " non-relativistic" 2D elec- 
tron gas (2DEG), in the systems of nodal fermions the 
single-particle density of states (DOS) can be strongly 
affected by disorder. Moreover, a greater variety of DOS 
behaviors and putative localization scenaria appear to be 
determined not only by the symmetry of disorder, but its 
strength as well^. 

In undoped graphene, the low-energy Dirac-like quasi- 
particle excitations reside near two conical points located 
at the momenta K — (47r/3V3a,0) and K' = — K (here 
a is the lattice spacing) in the inequivalent corners of the 
hexagonal Brillouin zone. These excitations can be con- 
veniently described in terms of the (retarded) bare Green 
function 



Gr{uj, p) = [(e -I- iO)% - Pf^jf,] 
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where fi 
1. The 



X, y and the fermion velocity is v — 
X 4 7-matrices %,x,y = (1 (8" cts,*! <8i 
(72,— j(73 (8" (Ti) act in the space of Dirac bi-spinors "0 — 
{^pK{^),'4'K{B),iJK'{A),i^K'{B)) composed of the val- 
ues of the K(K')-momentum components of the electron 
wave function on the A and B sublattices of the bipartite 
lattice of graphene. 

An important feature of the fermion spectrum de- 
scribed by Eq.(l) is its "chiral" (e — e) symmetry, 
which, if respected by disorder, facilitates a possible 
emergence of delocalized states at e = 0^. Indeed, un- 
like in the case of the conventional 2DEG prone to weak 
localization, the conductivity of graphene appears to be 
surprisingly robust, assuming values about three times 
higher than the ballistic result (cr = 2e'^ /-k^Tl) pertaining 
to the case of vanishing disorder^. 

In an attempt to explain this observation, a num- 
ber of proposals invoked the presence of an effective 
random magnetic field (RMF) causing a suppression of 
weak antilocalization that would otherwise have occurred 
due to an ordinary (chirally non-symmetric) potential 



disorderiSi^. In that regard, it should be noted that, be- 
ing a source of elastic scattering itself, a RMF, too, can 
be solely responsible for the onset of localization in the 
conventional 2DEG, the latter developing in accordance 
with the standard unitary scenario^. 

In this Letter, we will focus on the case of disorder 
which, in the framework of the Dirac equation, can be 
described as a random (and, generally speaking, non- 
Abelian in the valley subspace) vector potential with the 
Gaussian variance 



< i^(q) ® A^{-Ci) >= [1 ® lw{q) 



+ ((7+ (g) (T_ + ( 
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which is characterized by a parameter g = na'^ propor- 
tional to the areal density n of the RMF sources ("de- 
fects" ) with an internal length scale I. Beyond this length 
(q < 1//) the function w{q) decays algebraically 



(3) 



Although for any rj > the parameters g and I can be 
combined together, we choose to keep them separate. We 
also note that, owing to its effective nature, this RMF 
does not couple to physical fermion spin (if any) . 

With an eye on the problem of electronic transport in 
graphene, we will set out to explore the general case of 
vector disorder with a variable exponent 77 > 0. For 77 > 
and I ^ a the inter-node scattering does get suppressed 
{'w{l/l) ^ 'w{2K)), and the problem can be effectively 
treated as a single-node one (hence, in effect, Abelian). 
We will comment on the limitations of the single-node 
approximation below. 

The value of the exponent 77 is determined by the type 
of topological structural defects, such as discHnations 
(isolated pentagon- and heptagon-rings) with 77 = 1 and 
I ^ a, dislocations (pairs of adjacent pentagons and hep- 
tagons) with rj = and I given by the Burgers vector, 
etc.^. Moreover, a mathematically similar description 
was also proposed for the effect of thermal shape fluc- 
tuations of the graphene sheet ("ripples") whose origin 
can be traced back to the intrinsic thermodynamic insta- 
bility of 2D crystals. 

In recent derivations based on the theory of 2D elastic 
membranes, the effective gauge field was related to 
a fluctuating local height h of the graphene sheet with 
respect to the substrate (the authors of Refj^ derived 
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Af^ oc {dlh)"^, although a much stronger dependence, 
A-i_i oc {d^h)"^, can occur as welF). 

In what follows, we consider the more relevant version 
of RMF disorder described in Rcf7 where g ^ 1, I ^ 
bnm, and the exponent rj = 0.2 can be deduced from 
the previously studied asymptotic behavior of the height 
fluctuations of a 2D membrane of size L (< /i^ >cx L^+''). 

Despite the lack of immediate physical examples, we 
find it instructive to extend our discussion to the entire 
range of exponents < 77 < 2, for which the root mean 
square (RMS) value of the divergence-free part of the 
vector potential 



(2^ 



(4) 



{Ai oc hiL/l) increases with the system size, although 
the RMS of the physical magnetic field B = V x A 
and the corresponding energy density remain finite, Br/ oc 

It should be noted, though, that any value 77 7^ 1 
would be in conflict (see Ref.^) with the observed elec- 
tron density dependence of the conductivity of graphene, 
a ^ rie ^, if the ripples were to dominate over the other 
sources of elastic scattering. 

While the extensively studied Abelian and non-Abelian 
versions of the = problem are amenable to such pow- 
erful techniques as renormalization group, conformal al- 
gebra, and Liouville field theory, the case 77 > has, 
so far, only been studied for 77 = 1 and in the ballistic 
regime (e 3> l/^)^- In order to gain a preliminary in- 
sight, we first estimate the lowest-order correction to the 
(gauge-invariant) DOS 

Si^rjie) = Im I '^^^w{q){5^,i, - ^^)Tr[7^G7?X^^> P) 



1oGr{uj, p)^^jGr{uj, P + q)] 5 



1 - 77 e^-^'' 



which, incidentally, vanishes at 77 = 1, while at 77 = 
this correction {5i'o{t) ~ (jelne) is consistent with the 
previously obtained exact result 

u^{e) oc e^/^-i (6) 

where z = 1 + g for g < 2 (weak coupling) , whereas at 
strong couphng {g > 2) z = (85) — 1 — . 

It should be mentioned, however, that at 77 = and 
I ^ a the original model (2,3) becomes S'[/(2)-symmetric 
in the valley space {w{0) — 'w(2K)), in which case the 
power-law DOS features a universal exponent 1/7 — . 

The first order DOS correction (5) becomes compara- 
ble to the bare DOS (i/(°)(e) - e) at e - 5^/^''//, thus 
indicating the importance of higher-order corrections at 
still lower energies. Routinely, a further improvement 
would be sought out in the framework of the custom- 
ary self-consistent Born approximation (SCBA) for the 
(gauge-non-invariant) fermion self-energy 



Sfl(e,P) 



dq 



w{q) 



(2^)2Gfl(e,p-|-q)-i-SK(e,P + q) 
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It can be readily seen, however, that for 77 < 1 the SCBA 
solution remains finite at the nodal point, E/j(0, 0) — 
iT ~ ig^^^/l, whereas for 77 > 1 it is determined to the 

RMS of the RMF {T - AI,^^), both diverging with in- 
creasing L. A finite imaginary self-energy would then 
imply a finite DOS at e = (namely, f,,(0) ~ T\ InFaj ^), 
in a stark contrast with Eq.(6), which calls the validity 
of SCBA into question. 

It is well known that for 77 = the failure of SCBA 
stems from the fact that perturbative corrections repre- 
sented by diagrams with crossing disorder lines turn out 
to be as important as the non-crossed ones (which SCBA, 
in effect, accounts for)iS. However, for 77 > 1 the situa- 
tion appears to be more involved, since a singular (small 
angle) scattering tends to favor the so-called " maximally 
crossed" (including the nominally non-crossed) diagrams, 
as can be seen, e.g., from a direct comparison between 
the crossed and non-crossed second order corrections. 

Albeit not being a justification for SCBA, this observa- 
tion attests to a potentially singular behavior of the self- 
energy at small energies and/or momenta, which, in turn, 
would be indicative of a strong modification of the cor- 
responding quantum states by disorder that might even 
give rise to their complete localization. 

A further evidence to that effect can be gathered from 
the properties of the exact zero-energy states which can 
be explicitly contructed for any 77 > 0. In the trans- 
verse gauge representation of the random vector poten- 
tial Af^ = Ef^i^di/cj) describing a given RMF configuration, 
the exact (unnormalized) zero-energy states read 



7/)±(r) oc (l±7o) ( g_0(r) 



(8) 



Among the traditional objects of interest are such quan- 
tifiers of the (normalized) zero-energy wave functions' 
statisics as the inverse participation ratios (IPR) 



= < 



(/ e20(r')d2i./)n 



> 



(9) 



where the averaging over different RMF configurations is 
performed with the statistical weight 



=exp[- 



(27r)2 2w{q)^ 



(10) 



The IPR can be computed with the use of replicas, 
whereby one multiplies Eq.(9) by the denominator (nor- 
malization factor) raised to the power N — n, evalu- 
ates the resulting average, and eventually takes the limit 
— > 0, thus arriving at the result 

P„ = hm < / e2»^("-)+2^_^7^(r.)^2^ TT ^2 > 
•' i 

where for integer values of 77 the coefficient Crj{n) is a 
logarithmic function of L (e.g., Ci(7i) oc \nL/l). 
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The (nearly) quadratic (for all n) dependence on the 
system size is related to the groupping of the arguments 



(|r. 



I) which indicates a spontaneous symme- 



try breaking in the replica space. Furthermore, the ab- 
sence of rt-dependence (except, for, possibly, logarithmic 
factors) suggests that for rj > the zero-energy states 
tend to (pre)localize over a characteristic length of order 
I (possibly, up to a logarithmic factor). 

This result bears a certain resemblance to that found 
in the strong coupling regime of the previously studied 
77 = problems,. At weak coupling {g < 2) IPR ex- 
hibit a multifractal spectrum P„ ^ L^'^" governed by 
the anomalous dimensions A„ = {2 — gn){n — l) for g < 2 
and n < nc= (2/g)^/^ (A„ = 2n(l - l/ric)^ for n > ric), 
whereas for g > 2 this spectrum terminates (A„ — 
for all integer n), thus signaling the onset of a certain 
"freezing" transtion^^. 

Another signature of this peculiar (pre) localized be- 
havior can be gleaned from the correlation function of 
the (normalized) wave functions' amplitudes. Using the 
same replica trick, the latter can be cast in the form 



+i(dt,ip + dt,0)dt,(f> + -(cos2ip ~ cos2e)] (15) 
a 

The negative sign in front of the kinetic term of the field 
6 is a. salient feature of the bosonizcd ghost actioni^. 

Had all the integrals over the bosonic fields been con- 
vergent, one could have first integrated over 9, then 
shifted the field ip by icf) and, finally, integrated it out 
as well. The remaining integral over would then have 
all the operators ip in the integrand replaced by e"^, in 
agreement with Eqs.(ll,12). 

However, under a closer inspection one finds that the 
integral over appears to be controled by a modified 
weight 'P[(t)]ex-p{J d^qg^|(/)gp/(47r)^) and, therefore, di- 
verges for any r/ > (and g > 2 for rj — 0), thereby 
precluding one from the being able to readily integrate 
over ip and 9. 

Nevertheless, the problem can be circumvented if one 
first integrates over the disorder field, thus arriving at the 
" (sub)Sine-Gordon" action with a quadratic part whose 
diagonal terms remain positive definite for the modes 
with momenta q < g^/"^^ /I 



< 7/;2(ri)i/)2(r2) >=< 



g20(ri)g20(r2) 
(Jg20(r')rf2j.')5 



lim I 



27V-4 



exp[ 



-^(^-1)(^ 



r2| 



V 



> 



2771 



(2^ 



q\l+w{q))\^,\ 



(12) 



that clearly shows an unphysical ("negative norm") be- 
havior in the replica limit, which pattern is consistent 
with the replica symmetry breaking and a possible local- 
ization of the zero-energy states. 

From the formal standpoint, the origin of the above 
behavior lies in the divergence of the corresponding in- 
tegral over the disorder field for any rj > (at 77 = it 
only diverges for g > 2). In order to see that, we employ 
a supersymmetric formulation of the disorder averages 

exp(i J (fr[i)^^{id^ + A^) - e)7/> + 67^(^9^ + A^) - e)c]0l3) 

where X stands for an arbitrary product of local fermion 
operators and the path integral over bosonic ghost fields 
6 and c serves to the purpose of normalizing the fermion 
average by the partition function before averaging their 
ratio over disorder. 

A further progress can be achieved with the use of 
the 2D bosonization technique. Bosonization formulas^^ 
for the relevant bilinear operators {ijj'ip ~ a~^cos2(^, 
"0771 = iSfiudfj^ip and be ^ -a~^ cos26', 67^0 = i£^,i,d^9) 
involve two bosonic {ip and 9) and two auxiliary fermionic 
(^ and ^) fields, in terms of which the functional average 
(13) takes the form 

< X{i,, V.) > = j D[ip, 9, L i\VW\X{e''P)e-^ (14) 

and the bosonized action reads 

S\p, 9, 0] = / d'r[{d^^f - {d^9f + d^^d^i 



+ {w{q) - 1)\9,\' + \^,\' + 2w{q)^,9.,] (16) 

while its non- Gaussian part is given by the last term in 
Eq.(15). 

Albeit not being amenable to the renormalization 
group treatment for any 77 > 0, the resulting model can 
still be studied with the use of the variational method. 
In applications of this technique, the non-Gaussian part 
of the action is treated as a source of possible mass terms 
which cut off all the infrared divergencies. We choose a 
variational action in the form 



So[v^,9] - 



(2^ 



[q'{l + w{q))+ml)\^,\' 



+ iq\wiq)-l)+mi)\6 



2q^w{q)ipqt 



(17) 



where the mass parameters ?ti^ and me satisfy the cou- 
pled equations 



l_ 5^{S 
L2 ^ 5[ip,9Y 



So) 



>n=± 



D[ip, 9] cos2[ip,9]e~^" 
(18) 

For T] = the solution of Eqs.(18) reproduces the alge- 
braically vanishing DOS (6). The transition from weak 
to strong disorder and a concomitant (pre)localization of 
the zero-energy states (the states at e ^ are localized at 
arbitrarily weak disorder) which occurs at g = 2 can then 
be understood as a proliferation of unbound vortices of 
the field ip ii. 

At all < 77 < 2 Eqs.(18) allow for an (approximate) 
solution with the masses = = m?'{e) obeying 

the equation 



In 



777(e) 



d?ci q^w{q) 
(27r)2 [7772(e) +g2 



(19) 
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Solving Eq.(19) with a logarithmic accuracy, one obtains 
the mass 



1 



/|ln/e| 



(20) 



which reveals a chiral symmetry breaking and emergence 
of the chiral order parameter < tpip >oc< cos 293 > for 
e ^ 0. The solution (20) applies throughout the entire 
range of g and I, thus suggesting that ai rj > there is no 
counterpart of the freezing transition found at 77 
The corresponding DOS is given by the expression 



9.11 



1 



-Im < ipip >- 



1 



de ?7eP|ln/e|i+i/'; 



(21) 

which diverges at e — > but remains integrable for all 
rj. Interestingly, at ry = 1 the energy dependence (21) 
resembles that found in the unitary limit of random po- 
tential scattering^ (cf. Eq.(21) with the result reported 
in Ref.i3). 

It is worth noting that in the single-node approxima- 
tion the effect of (Abelian) potential scattering is not ex- 
pected to be localizing (which conclusion can be viewed 
as another manifestation of the Klein's paradox)^**. How- 
ever, the above results indicate that, by contrast, the ef- 
fect of a RMF on the Dirac fermions might be confining, 
consistent with the earlier observations^'"'^^. 

On the other hand, the authors of Ref.-lS argued that 
in the single-node approximation the effective description 
of undoped graphene can be achieved in terms of a pair 
of decoupled Wess-Zumino-Witten (WZW) or (if chiral 
symmetry is broken) non-linear a (NLcr) models, each of 
which has a topological term driving it to a stable fixed 
point of either symplectic (if the time-reversal invariance 
is intact) or unitary (as in the present case of a RMF) 
symmetry. 

However, the validity of the conclusions drawn in 
Refs.-^-fi hinges on the possibility of performing various 
disorder averages without encountering any divergencies 



in the underlying functional integrals. As follows from 
the previous discussion, such assumptions can only be 
justified in the weak-coupling regime {g < 2) of the short- 
range-correlated RMF (77 = 0), thus allowing for possible 
alternative scenarios for r] > 0. In that regard, it should 
also be mentioned that recent numerical simulationai^ 
seem to contradict the predictions of Ref.— pertaining 
to the non-chirally-symmetrical symplectic case. 

The above picture also suggests that, once the inter- 
node scattering has been put back in, the two WZW/NL(t 
models become coupled and the topological terms can- 
cel against each other, thereby restoring the standard 
localizing behavior, which in the RMF case would be de- 
scribed by the unitary NLcr-niodel (presumably, similar 
to that previously derived for non-relativistic spin- 1/2 
fermions with a gyromagnetic ratio equal two^^). 

Although the above scenario is believed to be rather 
generic, there might be such exclusions as the SU{2) 
( valley )-symmetric 77 = case of the model (2,3) where 
the zero-energy states remain critical (delocalized) for all 
g —. Moreover, it is also possible that a non-Gaussian 
{truly smooth) long-range-correlated RMF can support 
fermion trajectories which follow the lines of zero field 
and enable semiclassical percolation^^. This possibility 
lies well outside the domain of the Gaussian model (2,3), 
though. 

To conclude, in this work we applied a number of 
techniques to gain insight into the behavior of 2D Dirac 
fermions subject to a generic power-law correlated RMF. 
The available evidence suggests that the zero-energy 
states become localized, while the DOS diverges at e = 0. 

The results of this work can facilitate a better under- 
standing of the effects of ripples and topological struc- 
tural defects on electronic transport in graphene and 
other nodal fermion systems, such as the vortex line liq- 
uid phase of high- Tc cuprates. 
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